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Abstract 

We consider a family of translation-invariant quantum spin chains with nearest- 
neighbor interactions and derive necessary and sufficient conditions for these systems 
to be gapped in the thermodynamic limit. More precisely, let ^ be an arbitrary two- 
qubit state. We consider a chain of n qubits with open boundary conditions and 
Hamiltonian Hn{'ip) which is defined as the sum of rank-1 projectors onto ip applied to 
consecutive pairs of qubits. We show that the spectral gap of Hn{ip) is upper bounded 
by l/(n — 1) if the eigenvalues of a certain 2x2 matrix simply related to ip have equal 
non-zero absolute value. Otherwise, the spectral gap is lower bounded by a positive 
constant independent of n (depending only on ip). A key ingredient in the proof is a 
new operator inequality for the ground space projector which expresses a monotonicity 
under the partial trace. This monotonicity property appears to be very general and 
might be interesting in its own right. As an extension of our main result, we obtain 
a complete classification of gapped and gapless phases of frustration-free translation- 
invariant spin-1/2 chains with nearest-neighbor interactions. 
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1 Introduction 

Many properties of quantum spin chains depend crucially on whether the Hamiltonian is 
gapped or gapless in the thermodynamic limit. Ground states of gapped Hamiltonians 
are weakly entangled, as quantihed by the entanglement area law [1, 2, 3], and exhibit an 
exponential decay of correlation functions [4]. For such systems the ground energy and the 
ground state itself can be efficiently computed using algorithms based on Matrix Product 
States [5, 6, 7, 8]. On the other hand, ground states of gapless spin chains can exhibit 
drastic violations of the entanglement area law [9, 10, 11, 13], and computing the ground 
energy can be quantum-NP hard [14, 15]. Spin chain models studied in physics usually 
become gapless along quantum phase transition lines separating distinct gapped phases [16]. 
Deciding whether a given family of Hamiltonians is gapped or gapless in the thermodynamic 
limit is therefore a fundamental problem. 

In this paper we provide a complete solution of this problem for a class of translation- 
invariant chains of qubits with nearest-neighbor interactions. Let "0 G 0 be a fixed 
two-qubit state with ||-^|| = 1. Consider a chain of n qubits with open boundary conditions 
and define a Hamiltonian 

n—1 

HnW = Y, ll/^) (V’ki+l- (1) 

i=l 

Here each term is a rank-1 projector onto '0 applied to a consecutive pair of qubits. We 
shall refer to '0 as the forbidden state since the Hamiltonian penalizes adjacent qubits for 


2 










being in the state 'ip- As we will see in Section 2, the Hamiltonian Hni'ip) is frustration-free 
for any choice of ■0, that is, ground states of Hni'ip) are zero eigenvectors of each individual 
projector |'0)('0|j_i+i and the ground energy of Hni^) is zero. Futhermore, the ground state 
degeneracy of Hni^) is equal to n -|- 1 for almost all choices of xp. 

We are interested in the spectral gap separating the ground states and the excited states 
of Hni'ip) or, equivalently, the smallest non-zero eigenvalue of Hni'ip). To state our results, 
dehne a 2 x 2 matrix 


T,i, — 


( 2 ) 


(V>|0,1) (V^|1,1) 

-(V’|o,o) -(V^|i,o) 

Here |0), |1) is the standard basis of C^. As we will see, the matrix is crucial for un¬ 
derstanding the structure of the ground space of Hnpip). In this paper we prove that the 
eigenvalues of determine if Hn{'ip) is gapped or gapless^. Our main result is the following. 


Theorem 1. Letp) he an arbitrary two-qubit state. Suppose the eigenvalues ofT.^ have equal 
non-zero absolute value. Then the spectral gap of Hnpip) is at most l/(n — 1). Otherwise 
the spectral gap of Hni^ip) is lower bounded by a positive constant independent of n, which 
depends only on the forbidden state %p. 


We now motivate our choice of the model Eq. (1), highlight previous work on related 
models, and provide some intuition for why the eigenvalues of appear in the statement 
of the theorem. An informal sketch of the proof is provided in Section 1.1. Below we write 
'y{'ip,n) for the spectral gap of Hnigp). 

The family of Hamiltonians dehned in Eq. (1) includes some well-known quantum models 
as special cases. For example, choosing xp proportional to |0,1) — |1,0) (the singlet state) 
one can easily check that Hn{'ip) coincides with the ferromagnetic Heisenberg chain up to an 
overall energy shift. For this model Hnpip) has spectral gap 'y{'ip,n) = 1 — cos(7r/n) which 
decays as n~‘^ for large n [17]. Note that in this case T.iij is proportional to the identity matrix, 
so Theorem 1 gives an upper bound ''y{'ip,n) < l/{n — 1). Koma and Nachtergaele studied a 
one-parameter deformation of the Heisenberg chain known as the ferromagnetic XXZ chain 
with kink boundary conditions [17]. In this example ip is proportional to |0,1) — q'|l,0) for 
g > 0 and the spectral gap of Hn{ip) is given by 


■y{ip,n) = 1 — 2{qq ^cos(7r/n) 


for all n > 2, see [17] for details. One can check that is a diagonal matrix with eigenvalues 
Hi = {1 and H 2 = g(l + It follows that \hi\ ^ 1/^21 for any q ^ 1 and 

Theorem 1 asserts that Hniip) has a constant spectral gap. We note that in the two special 
cases considered above the Hamiltonian has a symmetry which enables an exact computation 
of the spectral gap. Such symmetries are not available for a general state ip. 

The exact results summarized above may suggest that the Hamiltonian Hniip) is gapless 
if -0 is a maximally entangled state and gapped otherwise. Theorem 1 demonstrates that 
this naive intuition is wrong. Indeed, choose ip proportional to a/ 1 — p|0, 0) -|- y/p|l,l) for 

^Although our definition of the matrix T.^ is basis-dependent, eigenvalues of 7/ are invariant under global 
SU{2) rotations. In particular, one can check that a transformation '0 —>■ (C/ 0 U)ip with a single-qubit 
unitary operator U maps T.^ to (dett/)“^ • UT^W. 
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some 0 < p < 1. Then the matrix has eigenvalues — p) and Theorem 1 implies 

that is gapless for all p as above. 

As a simple application of Theorem 1, we now map out the phase diagram of Hni'ip) 
restricted to the subset of real states -0 G 0 Using the Schmidt decomposition any 
real two-qubit state can be written as 


\i>±) = RiOi) ® R{e2) y/r^|o,o) ± ^|1,1) 


R{9) 


cos {9) sin {9) \ 

— sin (9) cos (6*) y ’ ^ 


for some 0 < p < 1/2 and 9i G [0,7r]. Since the spectrum of Hn{ip±) is invariant under a 
simultaneous rotation of all qubits, the spectral gap depends only on two parameters 02 — 0i 
and p. One can easily check that the eigenvalues of have equal non-zero magnitude iff 

4 

p>0 and sin^ (02 - 0i) < 2 + (^(i _ p))-i /2 • 

On the other hand, the eigenvalues of T^_ have equal non-zero magnitude iff either 
p = 1/2, or sin(02 — 0i) = 0 and 0 < p < 1. These conditions determine the gapless phase 
of the model for the special case of real states il). The gapped and gapless regions for -0+ as 
a function of p, 02 — 0i are shown in Fig. 1. A surprising feature is that the gapless phase 
occupies a hnite volume in the parameter space. In contrast, most of the models studied 
in physics only become gapless along phase transition lines which have zero measure in the 
parameter space. 



Figure 1: (Color Online) Phase diagram of the unfrustrated qubit chain where has 

real amplitudes. We use the parameterization \%Ij±) = R{9i)®R{92) [i/l — p|0, 0) ± ^^111,1)] 
and we show the gapped and gapless phases for -0+ as a function of 02 — 0i G [0,7r/2] and 
p G [0,1/2]. The phase diagram is symmetric under flipping the sign of 02 — 0i and under the 
transformation 02 — 0i —?• vr — (02 — 0i). The sector corresponding to V'- is not shown since 
it has a simple description: is in the gapless phase iff either p = 1/2, or sin (02 — 0i) = 0 

and 0 < p < 1. 

While it is possible to construct frustration-free translation-invariant Hamiltonians on 
qubits which are composed of projectors of rank 2 or 3, one can show that there are only a 
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handful of such examples. In the Appendix we describe them and for each we determine if 
the system is gapped or gapless. Taken together with our main result, this gives a complete 
classihcaction of gapped and gapless phases for frustration-free translation-invariant qubit 
chains with nearest-neighbor interactions. Note that the restriction to Hamiltonians which 
are sums of projectors is without loss of generality.^ 

There are several open questions related to our work. We do not know if the gapless 
phase of the model Eq. (1) can be connected to some known universality class of critical spin 
chains and what is the actual scaling of the spectral gap in the gapless phase. In particular, 
we do not expect that the the upper bound l/(n — 1) on the spectral gap in Theorem 1 
is tight. It is a challenging open problem to generalize our results to qudits, i.e., to map 
out the phase diagram of translation-invariant frustration-free spin chains for d-dimensional 
spins with d > 3. A natural analogue of the Hamiltonian dehned in Eq. (1) is 


n—1 

H„{m = ( 4 ) 

i=i 

where H is a rank-r projector acting on C'^ 0 C'^. It was shown by Movassagh et al [18] 
that such chains are frustration-free for any H and n whenever r < Their results 

also suggest that the Hamiltonian idn(n) may be generically frustrated for r > d^/4 and n 
sufficiently large. On the other hand, when r > d^/4, the Hamiltonian is frustration-free 
for certain special choices of H and n; examples include the famous AKLT model [19] (with 
d = 3, r = 5), the model based on Motzkin paths [11] (with d = r = 3), and the “Product 
Vacua with Boundary States” models [12] (with r = (d — l)(d -|- 2)/2). In general there is 
no efficient algorithm for testing whether idn (If) is frustration-free for a given n and there 
are indications that this problem may be computationally hard [20]. It is therefore natural 
to focus on the case r < d^/4, where the chain is guaranteed to be frustration-free. A next 
step could be to investigate the phase diagram of a chain of qutrits (d = 3) with projectors 
of rank r = 1, 2. 

Finally, if one moves from frustration-free one-dimensional chains to general two-dimensional 
systems, the problem of distinguishing between gapped and gapless phases of translation- 
invariant Hamiltonians becomes undecidable [21] which leaves no hope for mapping out the 
full phase diagram of such systems. 

1.1 Sketch of the proof 

Gapless phase. In Section 3 we consider the case when eigenvalues of T^j, have the same 
non-zero magnitude and prove that the spectral gap of is at most l/(n — 1). The 

proof uses a result of Knabe [22] relating the spectral gap of idn(V’) to that of the following 
Hamiltonian 

H°W = HM + \i^)mn,i (5) 

^Suppose instead we consider a frustration-free qubit chain H = where h has smallest 

eigenvalue zero (which can be arranged by adding a constant times the identity). Since H is frustration-free, 
it has the same null space as H' = X]r=i where 11 projects onto the range of h. Using this fact and 

the inequality cH' < H < ||/i||77', where c is the smallest non-zero eigenvalue of h, we see that H is gapped 
if and only if H' is. 
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which describes the chain with periodic boundary conditions. The other ingredient in the 
proof is a detailed understanding of the ground state degeneracy of We will see that 

is always frustration-free, but its ground state degeneracy can be smaller than that of 
In particular, if is not proportional to the identity operator then has a two- 

dimensional ground space whereas Hn{ip) has an n-|-1-dimensional ground space. Otherwise, 
I, then both Hamiltonians Hn{'il^) and have ground space degeneracy n -|- 1. 

We now sketch how these two ingredients can be used to prove the stated result. For ease 
of presentation we focus on the example considered above, where V’+ is of the form given in 
Eq. (3). Recall that the spectrum of Hn{iJ+) depends only on the two parameters 62 — 6 i,p. 
We can plot the ground state degeneracy of as a function of these two parameters. 

As described in the previous paragraph, this function takes the value 2 or n-|-1 depending on 
whether or not is proportional to the identity. The black lines in Figures 2(a) and 2(b) 
show the curves where the ground state degeneracy is equal to n -|- 1, for n = 10 and u = 50 
respectively. Everywhere else (0 < p < ^ and 62 — 61 G [0,7r/2]) the ground state degeneracy 
is 2. For reference we show the red and blue regions from Figure 1, which correspond to the 
gapped and gapless phases of the open boundary chain in the thermodynamic limit. As one 
might guess by looking at the Figure, the black curves become dense in the blue region when 
n —)■ 00 . If we consider a point 4^+ this blue region which does not sit directly on one of 
the black curves then the eigenvalue gap of Ff°('^+) is equal to its third smallest eigenvalue. 
However as n —?• 00 this point -0+ becomes arbitrarily close to a black curve, where the 
Hamiltonian has ground state degeneracy n -|- 1 and third smallest eigenvalue equal to zero. 
Using a bound on its derivative one can show that as a result the third eigenvalue of Ff°('0+) 
takes arbitrarily small values as u —)■ 00 . Finally, Knabe’s result implies that this can occur 
only if the spectral gap of Hn{'il>+) is at most l/(n — 1). This argument has to be modified 
slightly for states '0+ which, for some n, lie directly on one of the black curves and for general 
(complex) states -0. 

Gapped phase. In Sections 4,5 we prove that the spectral gap of line'll:) is lower bounded 
by a positive constant independent of n if the eigenvalues of have distinct magnitudes or if 
both eigenvalues are equal to zero. Our starting point is a general method for bounding the 
spectral gap of frustration-free spin chains due to Nachtergaele [23], see Lemma 3 in Section 5. 
To apply this method one has to manipulate expressions that involve the projector onto the 
ground space of Hn{'ip) which we denote The main technical difficulty that we had to 
overcome is a lack of an explicit expression for Gn which prevents us from straightforwardly 
applying Nachtergaele’s bound. Our proof is therefore indirect and is based on establishing 
some features of the ground space which allow us to control Gn sufficiently well. The key 
technical ingredient is a new operator inequality which expresses a monotonicity of the 
ground space projectors under the partial trace. More precisely, we show that 

Tr„(G'„) > Gn-i- (6) 

where the partial trace is taken over the u-th qubit. Using the fact that the Hamiltonians 
Ff„(0) are frustration-free one can easily check that Tr„(G„) and Gn-i have the same support, 
that is, Eq. (6) is equivalent to saying that all non-zero eigenvalues of Tr„(G„) are at least one. 
Our proof of this monotonicity property, presented in Section 4, applies to general frustration- 
free chains of qubits composed of rank-1 projectors. Neither translation-invariance nor the 
conditions of Theorem 1 are needed for the proof of Eq. (6). We note that Eq. (6) differs 
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(a) (b) 

Figure 2: (Color Online) Depiction of the ground state degeneracy of where '0+ 

is of the form given in Eq. (3). The black lines are curves in the {62 — Oi,p) plane where 
has ground state degeracy equal to n + 1. Here we plot the curves for (a) n = 10 
and (b) n = 50. For any point which does not lie on one of these curves (for 0 < p < ^ and 
6*2 — 6*1 G [0,7r/2]), the ground state degeneracy of is two. We also show the gapped 

(red) and gapless (blue) regions for the chain with open boundary conditions. As n —)■ cx) 
the black curves become dense in the blue region. 


from the well-known monotonicity property < Gn-i <8/. The latter follows trivially from 
the fact that Ffn('0) is frustration-free, whereas Eq. (6) holds for more subtle reasons. 

We proceed by showing that a quantum state which is completely mixed over the ground 
space of the n-qubit chain (i.e., proportional to the projector Gn) exhibits an exponential 
decay of correlations for certain local observables, see Lemma 4 in Section 5.1. In Section 5.2 
we use the decay of correlations and Eq. (6) to prove several “Region Exclusion” lemmas. 
Here we consider a partition of the chain into three or more regions and dehne local ground 
space projectors associated with each region. Loosely speaking, the Region Exclusion lem¬ 
mas state that the global ground space projector associated with the entire chain can be 
approximated by a certain operator built from the local ground space projectors. The latter 
are dehned on subsets of qubits where some of the chosen regions are excluded from the chain 
(hence the name of the lemmas). By repeatedly applying the Region Exclusion lemmas in 
Section 5.3 we arrive at the condition used in Nachtergaele’s bound, thus proving a constant 
lower bound on the gap. 


2 Structure of the ground space 

In this section we describe the ground spaces of and the Hamiltonians for the 

chain with open and periodic boundary conditions respectively (dehned in Eqs. (1,5)). 
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2.1 Open boundary conditions 

We first consider the Hamiltonian Hn{ip) for the chain with open bonndary conditions. We 
begin with the simple case where 0 V '2 is a prodnct state. It is always possible to 

choose the basis states |0) and |1) so that 


1 ^) = | 1 ) ® 1 ^^^) 


where 

|n) = c|0) + s|l), In-*-) = s*|0) — c*|l), and |cp + |sp = 1. 


For each i = 1,... ,n dehne an n-qnbit state \gi) = |0* *). Also dehne l^fo) 

For example, choosing n = 4 one gets 


l^o) 

\9i) 

\92) 

\93) 

1^4) 


V V V V ), 

V-^ V V V ), 

0 v-^ V V ), 

0 0 v-^ V ), 

0 0 0 ). 




Loosely speaking, the states Qi can be viewed as “domain walls” where |0) and |n) represent 
two different valnes of a magnetization. By direct inspection we see that go,...,gn are 
pairwise orthogonal gronnd states of Hni'il^). 


Proposition 1. Suppose s 7 ^ 0. Then the states go,... ,gn form an orthonormal basis for 
the ground space of Hni^). 

Proof. It snffices to show that the gronnd space of Hn{ip) has dimension at most n + 1. Dehne 
0 = 0 and 1 = n. Given any binary string x = (xi,..., Xn), dehne x = (xi,..., x„). Note 
that |0) and |1) are linearly independent since s 7 ^ 0. Therefore the states |x), x G {0,1}"' 
form a basis (non-orthonormal) for the Hilbert space of n qnbits. Snppose \g) is a gronnd 
state of Then l^f) = some complex coefficients a^- A simple calcnlation 

shows that 

i,i+l ("0 If?) "5 ^ ^ ) • • • ) 1) ^i+2 y • • • ■> ^n) 

x: {xi,Xi+i)={l,0) 

for any i = 1,... ,n — 1. On the other hand, i^i+i{'ip\g) = 0 since \g) is a gronnd state of 
Hnigp). This is possible only if a^, = 0 for all strings x that contain at least one consecntive 
pair (1, 0). Thns \g) belongs to a snbspace spanned by vectors where z = 0,..., n. 

This shows that the gronnd snbspace of Hn{'ip) has dimension at most n -|- 1. □ 

Now consider the case where is entangled. In this case we can still constrnct the gronnd 
space of althongh, in contrast with the prodnct state case, we are not able to obtain 

an orthonormal basis. The matrix T.^ dehned in Eq. (2) plays a crncial role. 

One can easily check that det (T^) 7 ^ 0 whenever is entangled and 


(-01 (/ 0 T^) = det (T^) (e 


( 7 ) 




where |e) = |0,1) — |1,0) is the antisymmetric state of two qubits. This shows that the 
ground space of hf2('0) = |'0)('0|i,2 is the image of the 2-qubit symmetric subspace under the 
map 1 ® T^. A similar characterization holds for Hn{'ip) with n > 2. In particular, dehne 

= ( 8 ) 

The following Proposition is a special case of a result presented in [24] (and has been used 
previously in, e.g., [25]). 

Proposition 2. Suppose det (T^) ^ 0. Then the ground space of is the image of the 

n-gubit symmetric subspace under the linear map 

Proof. Using Eq. (7) and the fact that M ® M\e) = det (M)|e) we get 

= |£){£|„+1 0 B, (9) 

for each j = 1,... ,n — l, where Bj is a positive operator acting on qubits in the set [n] \ {j, j + 
1}. From Eq. (9) we see that the nullspace of is equal to the symmetric 

subspace. The result follows since is invertible. □ 

Combining Propositions 1,2 and noting that the symmetric subspace of n qubits has 
dimension n -|-1, we conclude that the ground space of has dimension n -|-1 for almost 

any choice of (the only exception is when s = 0 and "0 is a symmetric product state). 

2.2 Periodic boundary conditions 

We now consider the Hamiltonian for the chain with periodic boundary conditions. It 

is well-known that is frustration-free for any choice of ip, see for instance [24, 26, 27]. 

However in this paper we will only need to deal with periodic boundary conditions in the 
case where is an entangled state. Accordingly, in this section we assume that det (T^) ^ 0. 
For any such we compute the dimension of the zero energy ground space of H°['ip). We 
will see that it takes different values depending on the choice of This contrasts with the 
open boundary chain which has ground space dimension n -|- 1 whenever ip is entangled. 
Here and throughout the paper we use the symbol ~ to mean proportional to. 

Proposition 3. Suppose Tfp ~ I. Then the ground space of Hfpip) has dimension n -|- 1. 
Otherwise, Hp^pip) has a two-fold degenerate ground space. 

Proof. Note that H'ppip) has the same rank as 

TpKWT? = t;'" (lo) 

where is given by Eq. (8). Both terms on the right-hand side are positive semidehnite 
and, by Proposition 2, the nullspace of the hrst term is the symmetric subspace. If ~ J 
then the second term in Eq. (10) can be written as |e)(e|„,i ® Bn where Bn is positive and 
|e) = |0,1) — |1, 0). Since this term annihilates every state in the symmetric subspace we see 
that in this case the nullspace of Eq. (10) is {n + l)-dimensional. 
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If is not proportional to the identity we show that there are exactly two states in the 
symmetric subspace which are annihilated by the second term in Eq. (10). We consider two 
cases depending on whether or not is defective (has only one eigenvector). 

First consider the case where has two linearly independent eigenvectors |ni), 1 ^ 2 )- Note 
that the last term in Eq. (10) projects qubits n, 1 onto a state 

w = (t;-*' 0 im ~ (t;' 0 /)|€) = (t;'|0>) |i> - (r;*|i)) |o>. 

The last equality makes it clear that \(f)) and |e) are linearly independent whenever is not 
proportional to the identity. Thus the nullspace of Eq. (10) is spanned by n-qubit symmetric 
states that are orthogonal to |0) on any pair of qubits. One can easily check that the only 
two-qubit symmetric states orthogonal to |0) are |ni ® vi) and \v 2 ® ^ 2 )- Likewise, one can 
check that the only n-qubit symmetric states orthogonal to | 0 ) on any pair of qubits are 
linear combinations of and \v 2 )®'^- Thus Eq. (10) has a two-dimensional nullspace and 

therefore the same is true for iL°('^). 

Next suppose is defective, i.e., has only one eigenvector. Let us work in a basis where 
| 0 ) is this eigenvector, so 



for some a, 6 G C with a 7 ^ 0. Then the last term in Eq. (10) projects onto a state 

\'t>) = {T;r» 0 im ~ (t;' 0/)|£> = b'\e) + a-|ii). 

Since a 7 ^ 0, the states |0) and |e) span the same subspace as |11) and |e). Therefore the 
nullspace of Eq. (10) is spanned by n-qubit symmetric states that are orthogonal to |11) on 
any pair of qubits. One can easily check that the only such states are linear combinations 
of |0)®"' and the n-qubit W-state 

| 100 ... 0 ) + | 010 ... 0 ) + ...+ 100 ... 01 ). 

Thus Eq. (10) has a two-dimensional nullspace and therefore the same is true for □ 

3 Gapless phase 

In this section we prove the hrst part of Theorem 1, namely. 

Gapless phase theorem. Suppose the eigenvalues of have the same non-zero absolute 
value. Then 7 (' 0 , n) < l/(n — 1) for all n > 2. 

Recall that 7 (' 0 ,n) denotes the smallest non-zero eigenvalue of iL„('0). In addition we write 
7 °(V’, n) for the smallest non-zero eigenvalue of the Hamiltonian with periodic bound¬ 

ary conditions, see Eq. (5). 

To prove the gapless phase theorem we use the following lemma, proven by Knabe [22], 
which relates the smallest non-zero eigenvalues of the chains with periodic and open boundary 
conditions. Knabe’s result, presented in Section 2 of reference [22], applies to more general 
frustration-free spin chains but here we specialize to the case at hand. 
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Lemma 1 (Knabe [22]). For all m>n > 2, 


1 




n — 1 
n — 2 




n — 



( 11 ) 


This lemma was originally proposed as a technique for proving that the periodic chain is 
gapped in the thermodynamic limit. This follows from the lemma if one can show that there 
exists a finite n for which the open chain has a gap strictly larger than l/(n — 1). Here we 
apply the lemma in the opposite direction. We use the following strategy which works for 
some (but not all) V’ satisfying the conditions of the gapless phase theorem. First we apply 
the argument sketched in Section 1.1 to show that 7°(V’, 'm) can take arbitrarily small values 
for large enough m. Then we apply Knabe’s lemma to infer that 'y{'ijj,n) < l/{n — 1) for 
any n > 2 since otherwise Eq. (11) would provide a constant lower bound on 7°('0,m) for 
all m > n, leading to a contradiction. Note also that 'y{'ijj,2) = 1 since = |'0)('^|. 

For some states '0 we are not able to use the above strategy directly; however in these 
cases we choose a state 0 which can be taken arbitrarily close to "0 for which the strategy 
can be applied. The result for -0 then follows by continuity. In order to handle these cases 
(and for other portions of the proof) we will need the following straightforward bound on 
how much the eigenvalues of Hni'ip) (or can change as V’ varies. Write 

ei{'ip,n) < e 2 {'ip,n) < ... < e 2 r>.{'ip,n) and el{'ip,n) < el{'ip,n) < ... < elrii'ipy’n) 
for the eigenvalues of ci-nd respectively. 

Proposition 4. Let 'ip and 0 satisfy 1101| = ||0|| = 1. Then 

\ej{^jJ,n) - ej{(j),n)\ < 2n\\^jJ - (j)\\ and |e°(0, n) - e°(0, n)| < 2n ||0 - 0|| 
for each j = 1 ,..., 2”. 

Proof. The proof of the two inequalities is almost identical so here we prove only the first one. 
We use the Weyl inequality for perturbed eigenvalues (see for example Corollary III. 2. 6 of 
reference [28]) which in this case says \ej{'ip,n) — 6^(0, n)| < ||if„(0) — if„(0)|| . To complete 
the proof we bound 


n—1 

WHni'lp) - Hn{(p)\\ < |||0)(0|i,i+l - |0)(0|i,i+i|| 

i=l 

= m ) - m m + (01) + (i^) + m m - (01)11 

<(n-1)11(10)-10)) ((01+ (01)11 
<2(n-l) 1110)-10)11 

where in the last line we used the fact that |||0) + |0)|| < 2 (since |0) and |0) are normalized). 

□ 


We now proceed to the proof of the gapless phase theorem. 
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Proof. First we claim that the eigenvalues of and absolute values of the eigen¬ 

values of are invariant under a transformation ip ^ {U ® U)^) where U is an arbi¬ 
trary single-qubit unitary operator. Indeed, let fj' = {U ® U)ijj. Then Eq. (1) implies 
= U^^Hn{'ip){U^)®^ and Eq. (7) implies = (det f/)"^ ■ UT^W. Thus the eigenval¬ 
ues of and magnitudes of the eigenvalues of T.^' do not depend on U. We shall use 

the freedom in choosing U to bring ip into a certain canonical form as dehned below. 

Proposition 5. For any {fj) G ® there exists a single-qubit unitary U such that 

{U®U)\ij) = (a + z/3)|0,l) + (a + z7)|l,0)+5|l,l) (12) 

for some real coefficients a, (3,^,6. 

Since the proof is rather straightforward, we shall postpone it until the end of this section. 
From now on we can assume that has the canonical form as in the right-hand side of 
Eq. (12). Substituting this canonical form into Eq. (2) one gets 

-(J-r-,))- ( 13 ) 

The eigenvalues of this matrix are a — i/3 and —{a — iy), with magnitudes and 

\/P 7 ^. By assumption of the theorem, the eigenvalues have the same magnitude and 
thus 7 = ±/3. We consider the two cases 7 = ±/3 separately and we show that 7('0, n) < 
in each case. 


Case 1: 7 = /? 

Fix n and let m > n be even. Setting /? = 7 in Eq. (13) and taking the square we see 
that T| ~ I. Since m is even we get T™ ~ I and therefore e'^l^ipi^) = 0 by Proposition 3. 
Let \cpm) be a normalized state which satishes 

\\(pm-7\\ < ^ 

and such that the eigenvalues of have different magnitudes and are both non-zero. 
This guarantees that is not proportional to the identity and det(T<^^) 7 ^ 0. Then, by 
Propositions 3 and 4, 

2 

7°(0m, m) = el{(pm, m) = el{cpm, m) - el{%p, m) < —. 


1 n — 2 

3 -77 ( 0 m, m) < 


2 

+ — 


n 


Applying Lemma 1 gives 

n — 1 n — 1 n — 1 m \n — 1 

for all even m> n, and using Propositions 3 and 4 again we have 

7(0,n) = e„+2(0m,’^) + (e„+2('0,n) - en+2{,(pm,n)) 

2n 


— f'n-|-2(0m, ^) T n 
n — 2 


< 


n 


m \n 


2n 

3 - 2- 


The result follows since m >n can be taken arbitrarily large. 
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Case 2: 7 = —/3 

Let a + il3 = where 9 and r are positive real numbers. Here r < 1 since V’ 

is normalized and r > 0 since the eigenvalues of are assumed to be non-zero. We hrst 
consider the case where 9 is irrational. In this case the convergents of the continued fraction 
expansion of 9 give sequences of positive integers {pj} and {qj} with 

^ - 9 < 4 . (19) 

gcd{pj,qj) = 1, and where {qj} diverges. Here we shall omit the hrst two convergents ob¬ 
tained by the standard continued fraction expansion, in order to guarantee that the sequence 
{qj} is strictly increasing and qj > 2 for all j. 

Dehne 0,- = —, and let 

3 qj ’ 


1^.) =re-'"(®i+f)|0,l)+re'"(®^+?)|l,0) +5|1,1). 


Then 


||Tj-Vi|| = rV 2 _ i| < r\f2\'K {9j - 9)\ < 




TT 


- 


(15) 


where we used the inequality |e®* — 1| < |(:c|, Eq. (14), and the fact that r < 1. Note that 

(5 


T*.- = 


0 ire *’^^9 


has eigenvalues Ei = and E 2 = . We have Ei ^ E 2 , which follows from the fact 

that 9j is not an integer, since qj > 2. Thus T^. is diagonalizable. Furthermore, E^^ = El\ 
and therefore (using the fact that it is diagonalizable) ~ /. Hence e°^{d!j,qj) = 0 by 
Proposition 3. On the other hand is not proportional to the identity since 9 is irrational, 
hence 7°('i4, qj) = 63 ( 41 , qj). Using these facts and Proposition 4 we have 


7°{^,Qj) = {el{^p,qj) - el{dlj,qj)) < 2qj ||Tj 


< 


2V2- 


IT 


Qj 


Now for all j such that qj > n we get 


, 1 n — 2 1 / n — 2\ 2y/27i 

7(4’, n) < - -^ - -J (4>, qj) < - 7 + 


n — 1 n — 1 


n — 1 


n — 1 


Qj 


and hence 7 ( 41 , n) < since the sequence {qj} diverges and the second term can be made 
arbitrarily small. 

It remains to consider the case where 9 is rational. In this case, for any e we may choose 
9' to be an irrational number satisfying 1^' — 6^1 < e. Letting 


) |0,1) + re'"(^'+t) |0,1) + (5|1,1) 
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we may now apply the above proof to get 7 ( 0 , n) = e„+2(0, n) < Now using Proposition 
4 we get 

7 (V’, n) = en+ 2 ( 0 , + (en+ 2 (V', - en+ 2 ( 0 , n)) < + 2n || |0) - 1-0) || 

<-h 2nrV27ie 

n — 1 

where in the second line we used the same reasoning as Eq. (15). Since e can be chosen 
arbitrarily small we get 7 (' 0 ,'n) < □ 

Finally, let us prove Proposition 5. 

Proof. Recall that a transformation -0 —)■([/ ® U)'if maps T.^ to (det U)~^ ■ UT.^U^. Here U 
is an arbitrary unitary operator. We shall choose a sequence of such transformations that 
bring T.^ into the canonical form dehned in Eq. (13) which is equivalent to Eq. (12). First, 
choose U such that |0) is an eigenvector of T.,p. This is always possible since any complex 
matrix has at least one eigenvector. Now we can assume that 


Tp = 


Pi 6 

0 P2 


for some complex coefficients pi,p25<5. Next choose U = where the phase 6 satishes 

Re(e*^(pi + P 2 )) = 0. This maps to and now we can assume that Re(pi + P 2 ) = 0. 
Finally, choosing U = diag(e*®, e“*®) one can map S to without changing pi,p 2 - Thus 
we can make 6 real. This brings into the canonical form dehned in Eq. (13). □ 


4 Monotonicity under the partial trace 

In this section we establish a relationship between the ground space projectors describing 
a chain of n and n — 1 qubits. Our result holds in a more general setting than considered 
elsewhere in this paper since we do not assume translation invariance. 

Let V’1,'02, • • • be an arbitrary sequence of normalized two-qubit states. For each 

n = 2,... ,m dehne a Hamiltonian 


n—1 

^n(V'l,'02, • • • ,- 0 ^- 1 ) = (16) 

which describes a chain of n qubits. This Hamiltonian is frustration-free for any choice of 
-01 ,... ,'0n-i (this follows directly from Proposition 6 given below). Let Qn be the ground 
space of iL„('0i,'02, ■ ■ ■ and Gn be the projector onto Qn- We adopt the convention 

Hi = 0 and Q\ = C^. 

First, we note that P ^n-i ® due to the fact that the considered Hamiltonians are 
frustration-free. This results in a trivial monotonicity property Gn < Gn-i 0 I- Below we 
prove that one also has a different type of monotonicity, namely Tr„(Gn) > Gn-i, where Tr„ 
represents the partial trace over the n-th qubit. 
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Lemma 2 (Monotonicity). For each n = 2,..., m one has 

> Gn-I. (17) 

Given the simplicity and generality of Eq. (17), one may be tempted to ask whether it 
holds for some trivial reason unrelated to the structure of the considered Hamiltonians. 
We have observed numerically that Eq. (17) can be false if G„_i and are chosen as 
projectors onto random linear subspaces Gn-i G and Gn G Gn-i ® C^, even if 

the dimensions of Gn-i and Gn match those of the ground subspaces of ■ ■ ■ 5 V'n- 2 ) 

and Hni'ipi ,..., V’n-i)- Thus any proof of the monotonicity property must exploit the special 
structure of the projectors Gn- In the absence of an explicit formula for Gn, one has to rely 
on some indirect arguments in order to derive Eq. (17). This partially explains why the 
proof of the lemma given below is rather cumbersome. 

Proof of Lemma 2. We use induction in n. The base of the induction is n = 2. In this 
case G 2 = I — |'0i)(V’i| and thus ^ 2 (^ 2 ) = 2/ — Tr 2 (|V’i)('?/’i|) > / = Gi- Here we used the 
fact that the partial trace of any two-qubit state is a density matrix which has eigenvalues 
at most one. We now prove the induction step. For brevity denote V’ = V'n-i such that the 
last term in Hni'Gi,-- - ,'ipn-i) is |V’)(^|n-i,n- 

First consider the case where ifj is unentangled, that is, ip = a ^ /3 for some single-qubit 
states a,G- In this case the result follows trivially without using the inductive hypothesis, 
since Gn-i ^ G Gn which implies Gn > Gn-i 0 and thus Tr(G„) > Gn-i- 

In the remainder of the proof we consider the case where ip is entangled (i.e., not a 
product state). Write the Schmidt decomposition of ip as 

\'G) = Vm\wo)\vo) + y/^\wi)\vi) (18) 

where {wpiUj) = (vpVj) = 6ij and po,Pi > 0 with Po + pi = 1. 

Let G^ = I — Gn- Obviously, Tr„(G„) = 2/ — Tr„(G;J;). Furthermore, the trivial 
monotonicity Gn ^ Gn-i ® implies that Tr„(G„) has all of its support on Gn-i, that is, 
TiniGn) = TTn{Gn)Gn-i = Gn-iTTn{Gn)- Dehue an operator 

Rn = Gn-lFrn{G:^)Gn-l- (19) 

The above implies that 

Tr„(Gn) = Gn-lFrn{Gn)Gn-l = 2Gn-l — Rn > (2 — || || )G„,-1. (20) 

Thus it suffices to prove that ||7?n|| < 1- 
Choose an arbitrary orthonormal basis 

9l, 92, - - - , Qr € Gn—2, {9oi\9ls) (21) 

Also choose an arbitrary orthonormal basis 

hi, ^ 2 , - - - ,hg G Gn—i, {hi\hj^ ^i,j- (22) 

In general the dimensions r and s of the spaces Gn -2 and Gn-i will depend on the states 
ipi,, ipn -2 but we will not need an explicit expression for them. We will however need to 
use the fact that s > r, which we now establish. The following proposition is a special case 
of the result presented in [18]; we include a proof here for completeness. 
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Proposition 6. Let be the dimension of Qn- Then > D^-i for all 2 < n < m. 

Proof. Recall our convention that Qi = C^, so Di = 2. On the other hand H 2 {'f>i) = 
|' 0 i)(V’i|i ,2 and D 2 = 3, which conhrms D 2 > Di. We now establish that Dn — -Dn-i ^ 
Dn-i—Dn -2 for all n > 3. This is sufficient to complete the proof since it implies Dn—Dn-i > 
D 2 — Di = 1 . 

Let 0 be a general state in Qn, with n > 3. Let 71 ,... 7 d„_i be an orthonormal basis for 
Qn-i and let ki, , kd „_2 be an orthonormal basis for Qn- 2 - We can write 

D„-i 

10) = Ao|7i)|0) + Ai|7i)|l) 

i=l 

for some complex coefficients The fact that the Hamiltonian Hni'Qi, ■ ■ ■ ,'Qn-i) is 

frustration-free implies 

Gn = {Gn-1 ® C^) n {Qn-2 ® Qn-l) 

and thus the dimension of Qn is the number of linearly independent solutions to the equations 


{kj 0 Qn-i\4>) = 0 for all j = 1 ,..., Dn- 2 - (23) 

This is a set of Dn -2 linear equations for the 2Dn-i variables {fi,z}- The number of linearly 
independent solutions satishes Dn > 2Dn-i — Dn- 2 , or equivalently Dn — Dn-i > Dn-i — 

Dn-2- n 

Dehne r x s matrices 

= {ga 0 wo\hi) and = {ga 0 wi\hi). (24) 

where Wo,Wi are dehned in Eq. (18). The trivial monotonicity Qn-i ^ Qn -2 ® implies 

M^Mo + MlMl = I,. (25) 


Here and below Ig denotes the identity matrix of dimension q. Furthermore, expressing 
Gn-i = \hi){hi\ one gets 

{gc\Tln-l{Gn-l)\gis) = {ga®WQ\Gn-l\gf}®WQ) + {ga®Wi\Gn-\\gf}®Wi) = {MqMI + MiM\)a,fi. 

Since Tr„_i(Gn-i) > Gn -2 by the induction hypothesis, we infer that 

MqMI + MiMl > Ir. (26) 

The usefulness of the matrices Mq, Mi comes from the following facts. Recall that we dehned 
Rn = Gn-lTrn{Gf^)Gn-l. 

Proposition 7. Suppose Q is entangled. Then the matrix of the operator Rn in the chosen 
basis {hi,..., hg} of Qn-i can be written as 

Rn = PoMI{pqMqMI + piMiM\) ^Mq + piM\{pqMqMI P piMiM\) ^Mi, (27) 

where Po,pi > 0 are defined by Eq. (18). 
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Proposition 8. Let r,s be arbitrary positive integers with s > r. Let Mq, Mi be arbitrary 
matrices of size r x s satisfying MqMo + M|Mi = Ig and MqMq +MiM| > Ir- Let po,pi > 0 
be any real positive numbers. Then the operator Rn defined by Eg. (27) satisfies ||-Rn|| < 1- 

Combining Propositions 7,8, the inequality s > r proved in Proposition 6, and Eqs. (25,26) 
one gets ||-Rn|| < 1- The lemma follows from Eq. (20). □ 

In the rest of this section we prove the above propositions. 

Proof of Proposition 8. Denoting x = Pi/po one can rewrite as 

Rn = mI^MqMI + xM^mIY^Mo + xMI{MoMI + xM^mIY^Mi. (28) 

By symmetry we can assume that x >1. Then Eq. (26) implies 

MoM(J + xMiMl = MqMI + MiM| + (x - 1)MiM| > R + {x - 1)MiMJ. 

Since the function fiy) = —1/y is operator monotone, we arrive at Rn < Sq + Si, where 

^0 = Mo^(/, + (x-1 )MiM/)-^Mo and = xM|(J, + (x - 1 )MiM/)-^Mi. (29) 


Hence it suffices to prove that Us'd + S'!]! < 1. From Eq. (25) one infers that ||Mo|| < 1 and 
||Tfi|| < 1. Since Mq has (s — r) fewer rows than columns, it must have at least this many 
linearly independent vectors in its nullspace. From Eq. (25) one infers that for any 0 G 
with II0II = 1 and Mo0 = 0 we have ||Mi0|| = 1. Thus Mi has at least (s — r) singular 
values equal to 1. Likewise, Mq has at least (s — r) singular values equal to 1. Note that this 
implies that {s — r) < r since MqMq is an r x r matrix with at least (s — r) eigenvalues equal 
to 1. Furthermore, conditions Eqs. (25,26) and the norm of Sq + Si are invariant under a 
transformation Mqj —)■ WMq^iV, where W and V are arbitrary unitary matrices. We can 
always choose W and V to bring Mi into a diagonal form such that the diagonal matrix 
elements of Mi are non-negative and non-increasing. Thus we can assume without loss of 
generality that 

Ml = [ D Orx(s-r) ] (30) 

where 


D = diag((ii, d 2 , ■ ■ ■ ,dr), 1 = di = ... = d(^s_r) > d(^s-r+i) > ■ ■ ■ > dr > 0. (31) 

(If s = r the above equation should read 1 > di > d 2 ■■■> dr > 0.). Here and below Oixq 
denotes an all-zeros matrix of size t x q. It follows that MiM\ = and thus 


Si = 


Rs—r) 
Orx (s—r) 


0 


(s—r)xr 


where Si = xD‘^{Ir + {x — 1)DR 


(32) 


and 


D diag(d(s_j.+l), d(^s—r+2)-! ■ ■ ■ i dr, 0ix(s—r))- 


(33) 


The above arguments also show that 0 G is in the nuillspace of Mq iff 0 has support on 
basis vectors i with di = 1. Thus we can assume wlog that 


Mq — [ Orx (s—r) M j , 


(34) 
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( 35 ) 


where M is some matrix of size r x r. Substituting Eqs. (30,34) into Eq. (25,26) yields 

M^M = Ir-D^ 


and 

MM^ >Ir-D^. (36) 

Using the polar decomposition of M and Eq. (35) we obtain the parameterization M = 
U{Ir — £1^)^/^, where U is unitary. Then Eq. (36) is equivalent to U{Ir — D‘^)U^ > h — 
or 

Ub'^U^ < £)2. (37) 

Using the dehnition of S'o one gets 


^0 = 


0(s—r) X (s—r) 
Orx (s—r) 


0 


(s-r)xr’ 

^0 


5 


where (J, + (x - l)Db~^M. 


(38) 


Combining Eqs. (32,38), it suffices to show that ||^o + >§i|| < 1- Using the chosen parame¬ 
terization of M one gets 


5o = {Ir - + {x- l)U^D^U)-\lr - ■ (39) 

Now Eq. (37) implies WD'^U > D^. Since f{y) = —1/y is an operator monotone function, 
it follows that 

{Ir + {x- l)U^D‘^U)-^ < {Ir + {x- l)b'^)~\ (40) 

that is, 

^0 + < {It - bb{Ir + {x- l)bb~^ + xb^Ir + {x - 1)1)^)-^ = I. (41) 

This proves that ||^o + >§i|| < 1 which implies ||S'o -|- S'!!! = 1 and thus ||i?n|| <1- □ 

Proof of Proposition 7. We first show that 

= A + B, where A = Gn -2 ® ip and B = ® C^. (42) 


Note that the two tensor products in A and in B refer to two different partitions of the chain. 
We use the following two general properties of the orthogonal complement (here W, V are 
subspaces of a Hilbert space) 

(>V + V)^ = >v^ n (43) 

(>V^V)^ =>V^^V + W®V^ + >V^®V^. (44) 


We have 

Gn = {Gn-1 ® C^) n {Gn-2 ® 

Applying Eqs. (43,44) gives 

Gn = {Gn-l ® C^)-^ -1- {Gn-2 8) 

= Gn-1 ® -f Gn-2 + Gn-2 ® b + Qn-2 ® V’ 

= A + B 
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where to get the last line we absorbed the second and fourth terms into the first, using the 
fact that {Gn -2 ® C 

Choose an arbitrary orthonormal basis 

ei,e2 ,...,eg e q = dim{g^_^) (45) 

From Eq. (42) one infers that is spanned by 

(ei,... ,eM) = (fi-i 0 <8) i/’) U (ei (g) no, ... ,eg (g) no) U (ei g) ni,..., (g) ni). 

where no,ni are the Schmidt vectors of tjj as dehned in Eq. (18). 

We now show that the Gram matrix F dehned by 


F 


p,g 


(CpICq) 


is invertible, which implies that ei,, cm are linearly independent. We note that F has the 
following simple form 


F 


Ir Bq Bi 

Bq Iq 0 

d 0 h. 


(46) 


where 

{Bg)a,i = {ga®g\ei®vg) = y/jh{ga®W:,\ei) 2) = 0,1. (47) 

Dehne B = \^ Bq i?i ], X = BB\ and Y = B^B. Note that X and Y have the same 
non-zero eigenvalues. Also note that F is invertible if none of these eigenvalues are equal to 
1 , since in this case 


(A-X)-i -B{h,-Y)-^ 

-B\lr-X)-^ ihq-Y)-^ 


(48) 


To show that F is invertible it therefore suffices to show that A — X is invertible. Using 
Eqs. (24,47) and the identity I = Gn-i + G^_i we get 


—BqBI + MqMI = Ir and —BM + MiM\ = h. 
Po Pi 


So 


Ir — X — Ir — BqBq — BiB\ — PqMqMq -|- PiMiM^. 


To prove that A — X is invertible we show that this operator is positive dehnite: 


PqMqM^ + PiMiMI > min (po,Pi)(AfoM(| -h MiM|) > min {po,pi)I > 0 


(49) 


where we used Eq. (26) and the fact that po,pi are both positive. This completes the proof 
that F is invertible and establishes that ei,... cm are linearly independent. 

Since we have shown that ei,... cm are a basis for we have 

M 

Gi= ^(r-‘)p.,|e,)(6,|. (50) 

p,q=l 
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Substituting Eqs. (50,48) into Rn = Gn-i^^n{Gn)Gn-i and noting that = ga®^^ with 
a = 1 ,... ,r are the only basis vectors of which are not orthogonal to Qn-i ® we 
arrive at 

r 

Rn=J2 ® Tr 2 |^)(^|)) G„_i. 

a,^=l 

Substituting Gn-i = \hi){hi\ and Tr 2 |'^)('^| = po|'n'o)('n'o| + Pi|'n'i)('?ni| into the above 
equation yields 

r 

{hi\Rn\hj) = '^{R- {po{h,i\ga ®wo) ■ {gis ® wo\hj) + pi{hi\ga ® wi) ■ {gp ®wi\hj)). 

q,/3=1 

Replacing the last two factors by matrix elements of Mq, Mi dehned in Eq. (24) and R — X 
by Eq. (49) one arrives at Eq. (27). □ 

5 Gapped phase 

In this section we prove the second part of Theorem 1, namely 

Gapped phase theorem. Suppose the eigenvalues ofT^ have different magnitudes or both 
eigenvalues are equal to zero. Then the spectral gap of Hn{R) is lower bounded by a positive 
constant independent ofn. 

Let us hrst consider the simple case when both eigenvalues of are equal to zero. Using 
the canonical form of R established in Proposition 5 and Eq. (13) one can check that this 
is possible only if \R) = {U ® U)|l, 1) for some single-qubit unitary operator U. Thus the 
Hamiltonian Hn{R) is a sum of pairwise commuting projectors and 'y{R,n) > 1 for all n 
which proves the desired lower bound. 

In the rest of this section we assume that the eigenvalues of have distinct magnitudes. 
In this case the eigenvectors of must be linearly independent. Let us hrst introduce some 
notation. Suppose S C [n] is a consecutive block of qubits. We shall write Gs for the 
projector onto the ground space of the truncated Hamiltonian 

{iji+ljCS 

obtained from Hn{R) by retaining only the terms fully contained in S. The projector Gs 
acts trivially on all qubits in the complement of S. Note that Gn = Gs in the case where S 
is the entire chain. 

Our starting point is a general lower bound on the gap of ID frustration-free Hamiltonians 
due to Nachtergaele [23]. Specializing Theorem 3 of Ref. [23] to our case one gets the following 
lemma. 

Lemma 3 (Nachtergaele [23]). Suppose there exists an integer r > 1 and a real number 
e < (r-|-l)“^/^ such that for all large enough n and for the partition [n] = ABG with \B\ = r 
and \G\ = 1 one has \\Gabc ~ GabGbcW < Then 

TiR, n) > (1 - e{r + 1)^/^)^ (51) 
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for all large enough n. 

Note that the right-hand side of Eq. (51) is a positive constant independent of n. Thus 
Lemma 3 reduces the problem of lower bounding the spectral gap of to that of upper 

bounding the quantity \\Gabc — GabGbcW- Our main technical result is an upper bound on 
this quantity that decays exponentially with the size of B. 

Theorem 2. Let /ii,/i 2 the eigenvalues of such that |/ii| < |/i 2 |- Define X = /i 2 /hi- 
Let c be the inner product between the normalized eigenvectors ofT^p. Consider any partition 
[n] = ABG such that \B\ = r. Then 

\\Gabc - GabGbcW < O (r'/ 2 |y|-./ 8 ) + ^ ^ (52) 


where the constant coefficients in 0{-) depend only on the forbidden state fj. ///ii = 0 then 
(52) holds with a formal replacement A = oo which sets the first term to zero. 

Note that |c| < 1 since the eigenvectors of T^j are linearly independent. Furthermore, 
since |A| > 1, the right-hand side of Eq. (52) is an exponentially decaying function of r. 
Therefore we can choose a constant r depending only on the forbidden state such that 

WGabc - GabGbcW < e = 2(r + l)i/2 

for all n > r. Substituting this into Lemma 3 one gets 


7(V’,n) > 


+ 1 ) 
4(r -|- 1) 


for all large enough n which proves the gapped phase theorem. 

In the rest of this section we prove Theorem 2. We shall hrst consider the case where 
is an entangled state (/ii ^ 0). The main technical difficulty we had to overcome is a lack 
of explicit formulas for the projectors Gabc^Gab, and Gbc- At a high level, our approach 
is to develop a set of identities relating the global ground space projector such as Gabc and 
the local ones such as Ga or Gab- These identities hold with a small error controlled by 
the size of the regions. Our proof of the theorem uses three identities of this type which 
are stated as “Region Exclusion” lemmas in Section 5.2. We use the Region Exclusion 
lemmas to decompose the operator Gabc ~ GabGbc in Eq. (52) into a sum of several terms 
and to show that the norm of each term is exponentially small in r, see Section 5.3. The 
proof of the Region Exclusion lemmas combines two ingredients: monotonicity of the ground 
space projectors under the partial trace (established in Section 4) and the fact that certain 
correlation functions in the ground space decay exponentially (established in Section 5.1). 

In Section 5.4 we consider the case where -0 is a product state (/xi = 0). In this case the 
orthonormal basis for the ground space constructed in Section 2.1 gives an explicit formula 
for the ground space projector. We use this formula to establish the Region Exclusion lemmas 
(for the /ii = 0 case) in a more direct way. The corresponding special case of the theorem 
then follows from the Region Exclusion identities. 

Before proceeding, let us establish some notation and conventions. Recall that a local 
unitary transformation 0 —)■ (f/ ® f/)^ preserves eigenvalues of Lfn(' 0 ) and maps to 
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(det U)-^-UT^U\ see Section 3. We shall choose the unitary U to hx one of the eigenvectors 
of T^. Specihcally, in the rest of this section we shall assume that 

T^|0)=/ii|0) and T^\v) = jj,2\v) (53) 

for some state 

|n) = c|0) + s|l), where |cp + = 1 and s > 0. 

Note that c is the inner product between the eigenvectors of T^, as dehned in the statement 
of Theorem 2. It will also be convenient to dehne a state 

= s|0) — c*|l). 

Given a set of qubits S and a single-qubit state \9) we shall write |6*)5 for the product state 
|5))®|S'I_ yyg shall write |6*)(6'|5 for the projector onto this state and |6*)(6'|;^ = Is — l^)(^|s- 

5.1 Correlation functions 

In this section we show that certain ground space correlation functions decay exponentially. 
Specihcally, dehne 


= Tr(G„|l)(l|i (g) |n-^)(n-^|j), 

(54) 

r(n) = Tr(Gn|n-^)(n-^|n) 

(55) 


where 1 < i < j < n. For notational convenience we have suppressed the dependence of 
these functions on the forbidden state ^jj. Our main result in this section is as follows. 

Lemma 4. The sequence {r(n )}„>2 is monotonically increasing and has a finite limit r(oo) 
such that 

0 < r(oo) — r(n) < O (nlAI”^"") for all n>2. (56) 

Furthermore, 

T{i-,j,n) < O ((j — i) ■ for all 1 < i < j < n. (57) 

Here the constant coefficients in O(-) depend only on the forbidden state fi. 

Proof. Let us dehne yet another correlation function 

a{i,j,n) = max (0| (|l)(l|i ® |n^)(n^|j) |0), (58) 

where the maximum is over normalized ground states, that is, ||0|| = 1. 

Proposition 9. 

cx(q j,n) < (59) 

1 — \c\ 

for 1 < i < j < n. 
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Proof. Define 


= |y^i) - (A*)nin^) 

Using Eq. (53) we see that (1|T^ = P 2 {M and (n'''|T^ = Therefore 

0 T^) ~ /i 2 (n-^l| - (A/ii)^(ln-^| ~ (n-^l| - (ln-^| ~ (e|, (60) 

where |e) = |0,1) — |1,0) (recall that ~ means proportional to). Comparing Eq. (7) and 
Eq. (60) one infers that is the forbidden state, that is, ip ^ 

Let 0 G be a normalized state (i.e., ||0|| = 1) for which the maximum in Eq. (58) is 
achieved, so 

= (01 (|l)(l|i ® |n^)(n^li) 10). 

Since 0 G by Proposition 2 it can be written \(p) = I ® ® T'^ ... ® where |x) 

belongs to the symmetric subspace. Using this fact and Eq. (60) we see that 

M(0^“i0)=O (61) 

for all integers 1 < i < j < n. Writing 

10 ) ~ |b)*l0o)N\* |n)i |0i) [n]\i 
and substituting this into Eq. (61) one gets 

s G(l|0o) - A^"*j(n^|01)) = 0 

which implies 

(011 (l^'^)(^'^lj) 101) = |^| 2 Q-q (001 (|l)(l|i) I0o)- (62) 

Using the fact that \(p) is normalized we have 

1 = (010) = (0*10*) + (01101) + 2Re (c(0* |0i)) . (63) 

We upper bound the magnitude of the third term using the Cauchy-Schwarz and the arith¬ 
metic/geometric mean inequality: 

|c(0*o|01)| < |c|-y(0(,|0(,)(0i|0i) < y ((0o|0o) + (01101)) • 

Substituting this into Eq. (63) yields 

1 > (1 - |c|) + wif,)) > (1 - |c|) 

and hence (0ol0o) < Using this fact and Eq. (62) we obtain 

and thus 

= (01 (|l)(l|i ® 10) = 5^011 (l^'^)(^'^lj) 101) < ^ 1^1 ■ 

□ 
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Now we are ready to prove Eq. (56). First, applying the Monotonicity Lemma (Lemma 2) 
to the left-right flipped chain yields Tri(G„) > Gn-i- Therefore 

r(n) = Tr(Tri(G„)|n-^)(n-^|„_i) > Tr(G„_i|n-^)(n-^|„_i) = r(n - 1), 

that is, T(n) is monotonically increasing. 

Inserting the identity decomposition I = |0)(0| -|- |1)(1| on the hrst qubit in Eq. (55) one 
gets 

r(n) = Tr(G„|l)(l|i 0 |n-^)(r;-^|„) Tr(G„|0)(0|i 0 |r;-^)(r;-^|„). (64) 

The hrst term in Eq. (64) is upper bounded by (n -|- l)(j(l,n, n) since we can decompose 
Gn = X]a=o \9a){9a\ usiug some orthonormal basis {ga} of Qn and use the fact that 

(fi'a|(|l)(l|i 8> \v-^){v-^\n)\9a) < a{l,n,n) 

for each individual state ga- The second term in Eq. (64) is upper bounded by r(n — 1), 
which follows from Gn < / 0 Gn-i- Thus 

r(n) < {n + l)o'(l, n, n) + T{n — 1). 

Proposition 9 implies a{l,n,n) = OdAI”^""), that is, 

0 < r(n) — T{n — 1) < 0(n|A|“^”'). 

This shows that r(n) has a hnite limit r(oo) at n —>■ oo. Summing up the series produces 
the desired bound Eq. (56). 

The proof of Eq. (57) follows a similar strategy. First consider the case i = 1, j = n. The 
same argument used above shows that 

r(l, n, n) < {n + l)cr(l, n, n) = (9(n|A|“^”). (65) 

Next suppose i >2 and j = n. Inserting the identity decomposition I = |0)(0| -|- |1)(1| 
on the hrst qubit, using the fact that < / 0 Gn-i, and noting that |l)(l|i < I, one gets 

T{i,n,n) < r(i - l,n - l,n - 1)-1-Tr(G'n|l)(l|i 0 |l)(l|i 0 |n-^)(n-^|n) 

< r(i — 1, n — 1, n — 1)-I-r(l, n, n). (66) 

This shows that 

n 

T{i,n,n) < T{l,k,k). (67) 

k=n—i-\-l 

Substituting Eq. (65) into this bound and summing up the series, we get 

r(f, n, n) < O ((n — i)|A|“^*'"'“*^) for all i = 1 ,..., n — 1. (68) 

(Here we included the case i = 1 which was handled in Eq. (65)). 

Finally, consider the case j < n — 1. Inserting the identity decomposition I = |n)(n| -|- 
|n“*“)(n“*“| on the n-th qubit, using the fact that Gn < Gn-i^I, and noting that < I, 

one gets 

r{ij,n) < T{i,j,n- 1) -f Tr(Gn|l)(l|i 0 0 |n-^)(n-^|n) < T{i,j,n- 1) -h r(i,n,n). 
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This shows that 

n 

T(h j,n) < (69) 

m=j 

Combining this with Eq. (67) leads to the desired bound Eq. (57). □ 

5.2 Region exclusion lemmas 

To perform manipulations with ground space projectors that involve several regions we now 
prove three region exclusion lemmas. These lemmas allow one to exclude one of the regions 
from certain operators built from ground space projectors. 

The hrst region exclusion lemma states that Gabc\v){v\bc ~ Gab <8) Ic\'v){v\bc- 

Region Exclusion Lemma 5. Let [n] = ABG with \B\ = j. Then 

\\{Gabc - Gab ® Ic) \v){v\Bcf < 0(|c|^) + 0(j|A|-^') 

Here the constant coefficients in O(-) depend only on the forbidden state xjj. 

Proof. Dehne P = Gabc and Q = Gab ® Ic- Using the fact that PQ = QP = P one can 
write the quantity we wish to bound as 

||(F - Q) |u)(u|sc 11^ < Tta{vbc\{P - Qf\vBc) = Pr{Q\v){v\Bc) - Pr{P\v){v\Bc)- (70) 
Define 

9{n,r) = Tr(G'n ■ ffi-r <8) |u)(u|®'’). (71) 

One can rewrite Eq. (70) as 

\\{P-Q) \v){v\Bc\f < 0{i + j,j) - e{i+ j + k,j + k), (72) 

where i = 1^41, j = \B\, and k = |C|. Representing |u)(u| = / — on the last qubit in 

Eq. (71) and using the monotonicity property Tr„(G„) > Gn-i from Lemma 2 one gets 

6{n, r) > 9{n — l,r — 1) — f{n,r — 1), (73) 

where 

^(n,r) = Tr(G„ • ffi-r-i ® ® |u-^)(u-^|). (74) 

From Eq. (73) we get 


j+k-l 

0{i+ j,j) <9{i+ j + k,j + k) + f{i + r + l,r) 

r=j 


and plugging into Eq. (72) gives 


j+k-l 

\\{P - Q)\v){v\Bcf < ^(^T’^ + If). (75) 

r=j 
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To complete the proof we now show that ^(n, r) has an upper bound which is exponen¬ 
tially small in r and does not depend on n. Partition the chain as [n] = A'B'B"C, where 
\C'\ = 1, \B'\ = [r/2j, \B''\ = l'r/2], and \ A\ = n — 1 — r. Using the fact that \v){v\b" < I 
we get 

^{n,r) < {vb'\Lb'\vb'), where Lb'= TTA'B''c'{Gn\v-^){v-^\c')- (76) 

Using the second part of Lemma 4 we have 

Tr{LB'\l){l\m) = T{'m,n,n) = 0{{n — for any m G 5'. (77) 

Note that n — m > r 12 for any m G B' . Let |0)(0|^/ = / — |0)(0 |b/. It follows that 

n—rl2 

lV(LB.|0){0|ij,) < ^ 1 V(Lb.| 1){1U) < 0(1) ^ (n - = 0(r|A|-'). (78) 

mdB' m=n—r 

Thus Lb' has almost all its weight on the basis vector |0 b/) and an exponentially small weight 
0(r|A|“'’) on the space orthogonal to |0 b'). Furthermore, the hrst part of Lemma 4 implies 
that Tt{Lb') = rin) = r(oo) — 0(n|A|“^"'). Combining this fact and Eq. (78) results in 

Lb'= t{oo)\0){0\b'+ E where ||i7|| < 0(r|A|“'’). (79) 

Therefore 

^(n,r) < {vb'\Lb'\vb') = r(cx))|(0|n)pl^'l0(r|A|"'’) < r(cx))|c|(^"^)-h 0(r|A|"^). (80) 

Finally, substituting this into Eq. (75) gives 

j+k—l oo 

\\{P - Q)\v){v\Bcf < ^(i + r-l-l,r) < ^^(z + r + l,r) < 0(|cp) + 0(j|A|"^). (81) 

r=j r=j 


□ 

The second region exclusion lemma concerns the operator Gabcd\v){v\c (recall that 
|n)(n|^ = I - |v)(v|c). 

Region Exclusion Lemma 6. Consider any partition [n] = ABCD with \B\ = j . Then 
II {Gabcd — | 0 )( 0 |^ ( 8 ) Gbcd) |' y )(' y|clP < 0 ( j | A | ^■^). 

Here the constant coefficient in O(-) depends only on the forbidden state xjj. 

Proof. For brevity denote P = Gabcd and Q = |0)(0|a ® Gbcd- Then 

||(F - Q)|t^)(r;|cir < Tr((P - Q)»(n|^). 

Taking into account that Gabcd Gbcd = Gabcd gives 

||(F - Q)|n)(n|cir < Tr(F|n)(n|^) + Tr(g|n)(nl^) - 2Tr(F|0)(0U ® |'y)(r^iy ■ 
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Substituting |0)(0 |a = I — |0)(0|^ in the last term results in 

||(F - Q)|u)(u|ciP < Tr(Q|u)(u|^) - Tr(P|u)('i;|^) + 2Tr(P|0)(0|ji ® |'y)('y|c)- 

Applying the Monotonicity Lemma (Lemma 2) one gets Tta{Gabcd) > Gbcd- This shows 
that 


Tr(Q|n)(n|c) - Tr(P|n)(n|c) = Tt{Gbcd\v){v\c) - Tt{Gabcd\v){v\c) < 0 
and therefore 

||(P - Q)|n)(n|cr < 2Tr(P|0)(0|:^ 0 |n)(n|p) < 2 ^ ^ Tr (Gabcd 11)(11m ® 

mSA m'£C 

One can recognize the last term as the correlation function T{m, m', n) defined in Section 5.1. 
Using the second part of Lemma 4 one gets 

OO 

ll(^’-Q)l«)(«lclt < E E O ((™' - < 0(l).5^ r(r-j)|A|-2’' = O(i|A|-^0. 

mSAm'sC r=j 

Here we denoted r = m' — m so that r > \B\ = j. We also used the fact that the number of 
pairs (m, m') with m & A and m' & G such that m' — m = r is at most r — \B\ = r — j. □ 

The third and hnal region exclusion lemma involves operators built from the ground 
space projectors as follows. Given any bipartition [n] = AB, where A and B are consecutive 
blocks of qubits, define 

Ga>b = Ga G) \v){v\b — Gab and Ga<b = |0)(0|a ® Gb — Gab- 
Region Exclusion Lemma 7. Consider any partition [n] = ABG with \B\ = j. Then 

\\Gab>c - |0)(0U ® Gb>c\\ < O + O (|cp/4) _ (§2) 

and 

\\Ga<bc - Ga<b ® |n)(n|c|| < O + O {\c\^/^) . (83) 

Here the constant coefficients in O(-) depend only on the forbidden state fj. 

Proof. We hrst show that the bound Eq. (83) follows from Eq. (82) and thus it suffices 
to prove the latter. To see this, consider horizontally flipping the chain so that the vertices 
previously labeled 1, 2,..., n are now n, n — 1 ,..., 1. The new forbidden state is = SWAP?/) 
where SWAP is the unitary transformation which interchanges the two qubits. The new 
matrix P^/ is proportional to From this we see that |0') = |n) and |n') = |0), and that 
A' = A. Using these facts we can see that Eq. (83) is just Eq. (82) applied to the left-right 
flipped chain. 

It remains to prove Eq. (82). Let 0 be any normalized state of ABG such that 

\\Gab>c ~ |0)(0U ® Gb>c\\ = \\{Gab>c ~ |0)(0U ® Gb>c)4>\\- (84) 
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Partition the region B as B = B'B", where \B'\ = [j/2j and \B"\ = \j/2]. Dehne states 

(j)~ = \v){v\c ■ (j), = \v){v\j^„ (S) \v){v\c ■ (j), = \v){v\b" ® \v){v\c ■ (j)- 

One can easily check that the above states are pairwise orthogonal, 

0 = 0- + 0-+ + 0++ and 1 = ||0||2 = ||0-||2 +||0-+||2_,_ ||0++||2_ 

We shall bound the contributions to the right-hand side of Eq. (84) coming from <p~,<p 
and separately. 

Let us start with (j)~. Using the dehnitions of Gab>c and Gb>c one gets 

Gab>c ' 4> =—Gabc'4> and |0)(0 |a ® Gb>c • 0 = —10)(0|^ ® • 0 • 

Since 0“ = |n)(n|^ • 0“ and ||0“|| < 1, this results in 

||(Gab>c — |0)(0|yi ® Gb>c)0 II < ||(Ga_bc — |0)(0|a <8 Gsc)|'i^)('i^lcl| < ^). (85) 

Here the last inequality follows from Lemma 6, where we set iA = 0. 

Next let us consider Using the dehnitions of Gab>c and Gb>c one gets 

Gab>c-4> ~^ = {Gab — Gabc)-4> ^ and |O)(O|a®Gs>c- 0 = |0)(0|yi®(G^ —G^c)'0 

It follows that 

||(Ga_b>c — |0)(0|yi ® Gb>c')0 "'"11 < ||(G^s/b» — |0)(0|a ® Gs/b")I'^)('^'Ib"II 

+ \\{Gab'B"c — |0)(0|yi ® Gb'B"c)\v){v\^„\\ (86) 

We can bound both terms in the right-hand side of the above equations using Lemma 6. One 
should choose the four regions in the statement of Lemma 6 as (H, B, G, D) = {A, B', B", 0) 
for the hrst term and (H, B, G, D) = (H, B', B", G) for the second term. This results in 

ll(G.4i.>c - |0){0U ® Gfl>c)r+ll < 0 (j‘G|a|GG) + 0 (j'G|a|-JG) = 0 (j'G|A|-""). (87) 

Finally, let us consider We have 

Gab>c'4>^~^ = {GAB — G abc) ' and |O)(O|^®Gs>c-0 ’'"’'' = |0)(0|ai® (G^ —Gsc) •0''"^- 
It follows that 

\\Gab>c ■ 0^’'’|| < \\{G(^ab')b"c — G(^ab')b")\'v){'v\b"c\\ < O -I- O (|cp/'^) . (88) 

Here we applied Lemma 5 choosing the regions in the statement of the lemma as (H, B, G) = 
[AB\B'\G). Likewise, 

II|0)(0|yi ® Gb>c'• 0''”'"|| < II(Gsc — Gs) • 0’''^|| <\\{Gb'b"c~Gb'b")\v){v\b"c\\ 

< O (j^/^|A|-^/^) + O (|cp/^) . (89) 


Here we applied Lemma 5 choosing the regions in the statement of the lemma as (H, B, G) = 
{B',B'',G). Substituting Eqs. (85,87,89) into Eq. (84) and using the triangle inequality one 
arrives at the desired bound Eq. (82). □ 


5.3 Proof of the gapped phase theorem 

Let us now prove Theorem 2 for the case where "0 is entangled, that is, /ii 7^ 0. In Section 5.4 
we will see how to modify this proof to handle the product state case /ii = 0. 

Partition region B as B = B'B”, where \B'\ = fy/2j, \B"\ = fy/2]. First we note that 

GabGbc — Gabc = {Gab — G abc)G bc (90) 

and 

Gab — Gabc = Gab ® + Gab < 8 ) |' y )(' y|c ~ Gabc = Gab>c + Gab <8 | n )( n |^. 

Here we used the notation from Lemma 7. Applying Lemma 7 to exclude region A from 
Gab>c one gets 

Gab — Gabc = |0)(0|a ® Gb>c + Gab ® + G (91) 

where denotes some operator such that 

||6.||<0(ri/2|A|-^/4)+0(|cr/"). 

Substituting the identity 

Gb>c = Gb 0 | n )( n|c — Gbc = {Gb — Gbc) — Gb ® b )(' y|c 

one gets 

Gab — Gabc = | 0 )( 0 |yi ® {Gb — Gbc) + {Gab — | 0 )( 0 fy 0 Gb) ® + G 

= |0)(0fy (g) {Gb — Gbc) ~ Ga<b'b" <8 + G- (92) 

Applying Lemma 7 to exclude region B" from Ga<b'B" one gets 

Gab ~ Gabc = |0)(0fy ® {Gb — Gbc) ~ Ga<b' ® \v){v\b" ® |'y)('y|c + G/2 + G- (93) 

Using {Gb — Gbc)Gbc = 0, ||Gyi<s'i| < 2, and Gbc = Gb"cGbc, one arrives at 

\\{Gab — Gabc)Gbc\\ G 2||(|n)(nfy" <8 \'v){'v\c)Gb"c\\ + \\g/ 2 \\ + IkrII- (94) 

Finally, partition B" = B'^B'^ with \B'l\ = [ \B''\ J and IH 2 I = ( \B''\ ] (so that each part 
has size k. r/A). Denote 

5 r /4 = \v){v\c{Gb"c — | 0 )( 0 fy " Gb’^c)- 

Applying Lemma 6 where the four regions are chosen as {A, B,G, D) = (5", i? 2 , C, 0) one 
concludes that 

II 5 ./ 4 II <0(r'/'|A|-^/^). (95) 

Replacing Gb"c in Eq- (94) by |0)(0fy" <8 Gb'^c and using Eq. (95) results in 

\\{Gab — Gabc)Gbc\\ G 2 || \v){v\b'' ■ | 0 )( 0 fy '' || + 2 || fy / 4 || + \\er/2\\ + fyr || 

< 2|(n|0)|l^"l + O (r^/'|A|-’'/®) + O (|cr/®) 

< +0(|cr/®) . (96) 

This completes the proof of Theorem 2 for the case when -0 is an entangled state. 
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5.4 Specializing to product states 

Finally consider the case /xi = 0. This implies that det(T^) = 0, that is, -0 is a product state. 
Using the notation from Section 2.1 write 

\^p) = \lv^), | n ) = c|0)+s|l), |n^) = s|0) - c*|l), 

where |cp + = 1. Here s ^ 0 which follows from the fact that /i 2 ^ 0. It is also easily 
checked that |0) and |n) are eigenvectors of corresponding to eigenvalues /xi = 0 ,p 2 = —s 
respectively. 

We now show that the region exclusion identities presented in Lemmas 5,6,7 (for the 
case of entangled 'ip) become exact equalities. Indeed, as was shown in Section 2.1, the 
ground space of has an orthonormal basis go, ■ ■ ■ ,gn, where l^fo) = and \gi) = 

| 0 *-iy-Lyn-i) ^ > 1 . Thus 

n 

Gn = J2\g^){g^\, (97) 

i=0 

We start with Lemma 6. Let [n] = ABCD be an arbitrary partition such that B and C 
are non-empty. We have to prove that 

(Gabcd — |0)(0|a ® Gbcd) ■ \v){v\q = 0. (98) 

Note that |n)(n|j ■ \gi) = \gi) for all i < j <n. This implies 

\gi){gi\ ■ |v)('y|c = \9i){9i\ ■ \'u){'u\c ■ \v){v\h = 0 all i e AB. 

Substituting Eq. (97) for Gabcd and using the above identity yields 

Gabcd ■ \9i){9i\ ' = |0)(0U ® Gbcd ■ 

ieCD 

since \gi) = |0)a ® \9i-\A\) bcd for all i G CD. This is equivalent to Eq. (98). 

Next consider Lemma 5. Let [n] = ABC be any partition such that B is non-empty. We 
have to prove that 

{Gabc — Gab 8) Ic) • \v){v\bc = 0. (99) 

Note that for any i E A one has 

\9i){9i\ = \9i){9i\A 0 \v){v\bc and \gi){gi\AB = \9i){9i\A ® 

On the other hand, \gi){gi\ ■ |v)(v|_bc = 0 for any i G BC. Therefore 

Gabc ■ \v){v\bc = \9i){9i\ ■ \v){v\bc = \9i){9i\A ® \v){v\bc 

jeAu{o} ieAu{o} 

Likewise, \gi){gi\AB ■ \v){v\bc = 0 for any i E B. Therefore 

{Gab<^ Ic) ■\v){v\bc = \9i){9i\AB ■\v){v\bc = \9i){9i\A<^\v){v\BC- 

ieAu{o} ieAu{o} 
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Comparing the last two identities one arrives at Eq. (99). 

Finally, note that the proof of Lemma 7 only uses Lemmas 5,6, and the fact that Eq. (83) 
is equivalent to Eq. (82) applied to the left-right flipped chain. In the proof of Lemma 7 
we use the fact that is invertible to establish this latter fact. In the case at hand is 
not invertible but since l"^) = we immediately see that Eq. (83) is just the left-right 

flipped version of Eq. (82). So both inequalities in Lemma 7 become exact equalities. 

The proof of Theorem 2 from Section 5.3 uses the Region Exclusion lemmas to establish 
the result. Since we have shown that each of these lemmas holds (with exact equality) for 
product states ip, we see that the proof of Theorem 2 also applies in this case if one formally 
sets A = oo in all error terms. 
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Appendix 

Qubit chains with higher rank projectors 

In the main body of the paper we considered qubit chains where n is rank-1 and i7„(n), 
defined in Eq. (4), is guaranteed to be frustration-free. Here we consider qubit chains where 
n is rank-2 or rank-3 and we determine which projectors H correspond to frustration-free 
chains. For each frustration-free chain we determine if the system is gapped or gapless. We 
say that i7n(n) is gapped if its spectral gap, denoted 7 (n, n), is lower bounded by a positive 
constant independent of n (otherwise it is gapless). We shall write Qn for the null space of 
Hnijl) and for the projector onto this space. 

The case where H is rank-3 is trivial, so we consider it first. In this case there is a unique 
two-qubit state |x) in the null space of H. If i7„(n) is frustration-free then there exists 
an n-qubit state -0 with reduced state on each pair of consecutive qubits i,i + 1 supported 
entirely on y. It follows that the rank-3 chain is frustration-free if and only if y = 6^ ® 6^ for 
some single-qubit state 6. Thus H = / — \6){6\®‘^ and i7„(n) is a sum of pairwise commuting 
projectors. This shows that Hnijl) has unique ground state |6')®"', and its eigenvalue gap is 
equal to 1, for all n > 2. 

The rank-2 case is slightly more interesting. There is a trivial case where H is a 1-local 
projector, i.e., n = P®Jorn = /®P; in that case 7 (n, n) = 1 for all n >2. The following 
theorem handles all other cases. 

Theorem 3. Suppose H is a two-qubit, rank-2 projector which cannot he written as I ® P 
or P ® I for some projector P. Then the dimension of the null space of Hn JA) satisfies 
dim(^„) = dim(^ 4 ) G {0,1, 2} for all n > 4. Moreover, exactly one of the following holds: 

1 . Q 4 , = span{|aaaa)} for some single-qubit state a. 

2. Q 4 = span{|aaaa), \(3(3(3(3)} for some linearly independent single-qubit states a,[3. 

3. Q 4 = spaxi{\a(3a(3), \(3a(3a)} for some linearly independent single-qubit states a,j3. 

4 - G 4 = span{|aaaa), |a-*-Q!a!a) -7 f\aa-^aa) -|- p\aaa^a) + p\aaaa^)} for some or¬ 
thonormal single-qubit states a,a-^ and non-zero / G C. 

5. Q 4 is empty. 

In cases 1, 2, 3, and 4, the Hamiltonian Hnjl) is frustration-free for all n > 2, whereas 
in case 5 it is frustrated for n > A. P„(n) is gapped in cases 1, 2, and 3, and it is gapped 
in case 4 if \ f\ //' I/I = 1 case 4 then the Hamiltonian is gapless, with spectral gap 
upper bounded as 7(n,n) < (1 — cos(7r/n)). 
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To establish lower bounds on spectral gaps we shall use Nachtergaele’s criterion [23]. 
Recall that Lemma 3 states this criterion for the case where If is rank-1. More generally 
if iLn(n) is frustration-free for n > 2 then the same bound holds for its spectral gap (with 
7('0, n) and 7('0, r -f 1) replaced by 7 ( 11 , n) and 7 ( 11 , r -|- 1) in the statement of the Lemma). 
For our purposes it will be sufficient to use the following weaker version of the bound. 

Lemma 5. Suppose ifn(n) is frustration-free for n > 2. Let [n] = ABC with \C\ = 1, 
\B\ = r, and |y4| = n — r — 1. Suppose there exist constants 0 < 5 < 1 and K > 0 such that 
for all sufficiently large n we have \\Gabc ~ GabGbcW < K5'^. Then Hn(jA) is gapped. 

Proof. We can always choose r so that K6^ < Plugging this choice into Nachtergaele’s 

bound with e = obtain 7 ( 11 , n) > 7 ( 11 , r -|- l)(4r -|- 4)“^. □ 

Proof of Theorem 3. We hrst establish that the range of If is spanned by two linearly 
independent states 0, "0 which are both entangled. It is easy to check that the only two 
dimensional subspaces of ® which contain only product states are of the form 
or ® X for some single qubit state y. By the hypothesis of the theorem If cannot be 
written as / ® P or P ® J, which implies range(n) does not have this form; thus it contains 
at least one entangled state 0. Let v G range (If) be linearly independent from <p. It is easy 
to see that we can always choose = cf zu for some nonzero z G C so that 0 is entangled. 

So range(n) = span{0,0} where 0,0 are both entangled (equivalently, and are 
both invertible). Furthermore, an n-qubit state y is in the null space of P„(n) if and only 
if it is in the null space of | 0 )( 00 j+i and | 0 )( 00 i+i for alH = 1 ,..., n — 1 . 

To complete the proof we now suppose that ^4 is nonempty and we consider two cases 
depending on whether or not Tff^T^p has two linearly independent eigenvectors. The theorem 
follows directly from the following propositions which handle the two cases. 

Proposition 10. Suppose A is a two-gubit projector such that range(n) = span{0,0} where 
are both entangled. Suppose thatTff^T.^j; has linearly independent eigenvectors {a,0} and 
that Qi is nonempty. Then one of the cases 1,2, or 3 from Theorem 3 occurs. Moreover, 
P„(n) is frustration-free for all n >2 and it is gapped. Its ground space dimension satisfies 
dim(^„) = dim(^ 4 ) for all n > 4. 

Proposition 11. Suppose A is a two-gubit projector such that range(n) = span{0,0} where 
0, 0 are both entangled. Suppose that Tf^T.^ has only one linearly independent eigenvector 
{a} and that ^4 is nonempty. Then case 4 from Theorem 3 occurs. Moreover, Hn{A) is 
frustration-free and has ground space dimension egual to 2 for all n >2. It is gapped if and 
only if I/I 0 1; if\f\ = 1 then j{A,n) < (1 - cos( 7 r/n)). 

□ 


In the remainder of this section we prove Propositions 10 and 11. 

Proof of Proposition 10 We hrst establish that Q 2 is spanned by 1 ® T^|a)|a) and 1 ® 
T^|0)|0). These states are linearly independent (since a,j3 are). Since Q 2 is 2-dimensional 
it suffices to establish that it contains both of these states. Clearly each is orthogonal to 
|0)(0| since |0) ~ (e|/ ®Tjf^ (recall |e) = |0,1) — |1,0)). To check that they are orthogonal 
to |0)(0| use the fact that (0| ~ {e|/ ® and that a,0 are eigenvectors of 
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Now consider n = 3. Since is nonempty, there exists a state |x) G Qs, which by 
Proposition 2 can be written as |x) = 1 ® ® T||s) for some state |s) in the 3-qnbit 

symmetric snbspace. Since the chain is frnstration free, the hrst two qnbits of y have all of 
their snpport in Q 2 - This implies |s) = a|a)®^ + b\f3)^^ where a,b are not both zero. By 
symmetry we assnme withont loss of generality that a 7 ^ 0. Next, imposing orthogonality to 

I'#’) {'#’12,3 gives 

a\a) ((£|r* ® a)) + 6|/3) ® T^^Tl\l3,l3)) = 0. 

Since a,/? are linearly independent, both terms mnst be zero. Using the fact that any 
two-qnbit state orthogonal to e is symmetric, and that a 7 ^ 0 , we get: 

• T^\a) is an eigenvector of and 

• If 6 7 ^ 0 then T^\(3) is an eigenvector of 

Now recall that a,/? are linearly independent eigenvectors of Note that 

not proportional to the identity (since this wonld imply that 0 is proportional to 0 ) and 
therefore a,l3 are the only eigenvectors of Hence T^\a) is proportional to one of the 

states q;,/ 9, and if 6 7 ^ 0 then the same holds for T^\(3). However, since is invertible, it 
cannot be the case that T^|a) ~ T^\(3)- Pntting this together we see there are 4 snbcases to 
consider (below we show that the hrst three correspond to cases 1., 2., 3., from the statement 
of the theorem and that the fonrth does not occnr): 

Case (i): T^\a) ~ |a) and T^\(3) = c\a) + d\l3) where c, d are both non-zero. 

In this case T^|/3) is not an eigenvector of which implies (by the facts established 

above) that 6 = 0 and thus |x) ~ |a)®^ is the only state in Likewise, la)®” is the only 
state in Qn for n > 3 (since it is the unique n-qubit state such that any three consecutive 
qubits i,i + l,i + 2 have all of their support on ^ 3 .) This establishes that we are in case 
1 from Theorem 3. Note that for any partition [n] = ABC with \B\ > 1 we have Gabc = 
|a)(a|®” = GabGbc and so the conditions of Lemma 5 are satished (with 6 = 0) and Hn(Jl) 
is gapped. 

Case (ii): T^|a) ~ |a) and T^\(3) ~ |/d) 

In this case (by the facts established above) Q 2 is spanned by |a)|a) and |/3)|/d). This implies 
that Qn is spanned by |a)®” and |/d)®” for all n > 2, so we are in case 2 of Theorem 3. One 
can easily construct an orthonormal basis and confirm that ||G„ — |a)(a|®” — |/5)(/d|®”|| = 
0(1 (a 1/3)I”). Using this expression three times and the triangle inequality we get \\GabGbc~ 
GabcW — -^l(o:|/3)|'^' where K is a constant, for any partition [n] = ABC. Hence the 
conditions of Lemma 5 are satished (with 5 = |(a|,3)|) and Hn(jA) is gapped. 

Case (iii): T^|a) ~ |/3) and T^|/3) ~ |a) 

In this case Q 2 = span{|a)|/3), |,3)|a)}, which implies that for all n > 2, Qn is spanned by 
product states |z/i) = \a)\(5)\a) ... |/3) and \i' 2 ) = |,3)|a)|,3)... |a) (the last tensor product 
factors are instead |a), |/3) respectively if n is odd). This shows that we are in case 3 of 
Theorem 3. Now constructing an orthonormal basis we see that ||G„ — \vi){yi\ — |z/ 2 )(z/ 2 ||| = 
0(|(a|/3)|”). Letting [n] = ABC and using this expression we get \\GabGbc ~ GabcW — 
7L|(a|/S)|l'®l for some constant K. Hence the conditions of Lemma 5 are satished (with 
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6 = |(a|/3)|) and Hn(Il) is gapped. 

Case (iv): T^\a) ~ 1/3) and T^\f3) = c\a) + d\l3) where c,d are both non-zero. 

In this final case 6 = 0 and the only state in is |x) = 1 ® ® T||a)|a)|a) ~|a)|/3)T^|/3). 

Any non-zero state k G ^4 must have its first three and last three qubits in the state y. This 
implies in particular that k is a product state. Furthermore, the second qubit of k must be 
in the state (3 (look at the first three qubits) and also in the state a (look at the last three 
qubits), which is impossible. Therefore Qi is empty, which is a contradiction. So case (iv) 
does not occur. □ 

Proof of Proposition 11 Define an orthonormal basis {|(i) = |a), |i) = lo"'')}. Since 
| 0 ) is the only eigenvector of in this basis we have 

= ( / ^ ) ( 100 ) 

for some c,d^C with c 7 ^ 0 (since det 7 ^ 0 ). 

First consider n = 2. We claim that 1(8)T^|6)|6) and 1(8)T^ (|0)|i) + |i)|0)) span ^ 2 - To 
see this note that these states are linearly independent and orthogonal to To show 

that they are also orthogonal to |0)(</>|, use the fact that (01 ~ (e|l (g) and Eq. (100). 

We now show that our assumption that ^4 (and therefore also Q 3 ) is nonempty implies 
1 ® T,/i ® r||0,0,0) is in Q^. To reach a contradiction, assume 1 (g) (g) T^|0,0,0) 0 Q 3 . 

By Proposition 2 any state OJ E Q 3 satisfies |a;) = / ® (g) T||s) for some |s) in the 

three-qubit symmetric subspace. The first two qubits of |s) must be supported entirely in 
^ 2 - Using the form of Q 2 derived above, we see that this implies |s) is a superposition of 
the symmetric Hamming weight zero and one states (with respect to the 0,1 basis). Since 
l(g)T^(g)T||0, 0, 0) 0 Qs (and since Q 3 is nonempty), we have shown that Qs is one-dimensional 
and contains a single state u of the form 

\uj) = l(^T^(^ T| (a| 0 , 0 , 0 ) + 6 |i, 0 , 0 ) + 6 | 0 , 1 , 0 ) + 6 | 0 , 0 , i)) 

where 6 7 ^ 0. By assumption there exists a state in Q 4 , which must have its last three and 
first three qubits each supported on Q 3 , i.e., it must be of the form |di)|ci;) = \uj)\ 62 ) for some 
single-qubit states 61 , 62 - However this is impossible since u is not a product state whenever 
6 7 ^ 0. Having reached a contradiction we conclude T||0,0, 0) is in Q 3 . 

Now 

2 , 3(011 T||0, 0,0) ~ |0) ((e|T^ 0 0)) = 0 

implies that Tp|0) is an eigenvector of By assumption this operator has only one 

eigenvector 0, so T|/,|0) ~ |0). Thus |0) is an eigenvector of both and T^. Therefore, in 
the 0,1 basis we have 

rp f X y \ f q r \ 4 f sx sy - rz \ 

^*=[0 z ) ^*^[0 s ) 0 j 

for some x, y, z, q,r,s E C. Comparing with Eq. (100) we see that sx = qz ^ 0. Now 

| 0 )=a; 06 ,i)-^ 0 i, 6 ) + 2/0i,i) (101) 
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and 


10) =g*|0,i) -s*|i,0) +r*|i,i) = {q*/x*) (a:*|0,i) -; 2 *|i, 0 )) +r*|i,i). (102) 

Comparing Eqs. (101),(102) we see that range(n) = span{|0), 10)} = span{|l, 1), |z/)} where 

/7-rTm (l”-t-/lt0» (103) 

vi + l/r 

and / = z*/x* (note that x,z are both nonzero since sx = qz Moreover these states 

are orthonormal so 11 = |i, 1)(1,1| + \v){y\. 

By Proposition 2, a basis for the zero energy gronnd space of X]* is given by 


Tf Y. 1^)’ J=0...n 

2:e{0,l}®",|2|=i 


where T^, ~ diag(l, /) and = 1 ® Ty ® ® ... ® T^~^. It is then easy to see that the 

only states in this space orthogonal to Yhi |1) l)(l 5 ^ki+i fh® basis vectors corresponding 
to j = 0,1 and linear combinations thereof. These two states span Qn for all n > 2; this 
shows that we are in case 4 of Theorem 3. 

We now show that the system is gapped if |/| 0 1. Withont loss of generality we may 
assnme |/| < 1; note that if |/| > 1 we may relabel the qnbits 1, 2, 3,..., n as n, n — 1,..., 1 
(flipping the chain left to right) which sends / —)• f~^. Define = |0)(0|®” + |x^’"')(x^’"'l 
where 


W® = VWTF 



i—1 


16 ^ 




1 0® 


1 < A; < n. 


(104) 


Noting that the gronnd space projector is = |0)(0|®"' + ||x"'’’^ll and that 


Ix’‘'”)(x"'”l - lx‘'”)(x' 


k.n I 


< 2 


lx”'") - lx‘'”)|| < 



1/2 


2 (i-l) 


= 2i/r 


we obtain 

\\G^n - GnW < 1 - + |||x"’")(x"’1 - Ix"’")(x"’”||| < l/k" + 2|/k < 3|/|k 

Now let a partition [n] = ABC be given with l^l > \B\ and \C\ = 1. Using the above 
bonnd three times, the triangle ineqnality, and the facts that HGjill < 1 and < |/|'^', 
we obtain 

WGabGbc - GabcW < WG'abG'bI: “ GabcII + 9|/I'®'- (105) 

Using explicit expressions for and we see that G^^^G^^q = and the 

first term above is zero. This shows that the conditions of Lemma 5 are satisfied with 6 = |/|, 
and iLn(n) is gapped. 

Finally, we establish that iLn(n) is gapless if |/| = 1. Using Eq. (103) we see that 
iL„(n) commntes with the total Hamming weight operator X]r=i lu(l|* therefore block 

diagonal with a block for each Hamming weight 0,..., n. The spectral gap 7 (n, n) is npper 
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bounded by the smallest nonzero eigenvalue within any given block. The matrix of the block 
with Hamming weight 1, in the orthonormal basis {\ei) = 1 0"'“*) : i = 1,... ,n}, 

is given by 

1, iii=j = loii=j = n 

2, if 2 < i = j < n — 1, 

-1, if \i- j\ = 1 
0, otherwise. 

This is 1/2 times the Laplacian matrix of the path graph of length n, and its spectrum is 
known. In particular, its smallest nonzero eigenvalue is (1 — cos(7r/n)), which upper bounds 
7(n,n) if I/I = 1. □ 
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